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ABSTRACT 

Previously published mtDNA sequence data have suggested 
that an undescribed species of Physa (“Species A”) may inhabit 
the swamps and ditches in the southeastern Atlantic coastal 
plain. These snails are characterized by slender shells and dark 
bodies, but are otherwise similar to the more widely 
distributed P. pomilia. Mate choice tests revealed significant 
sexual isolation between Species A and P. pomilia , with homo¬ 
gametic pairings of P. pomilia five times more frequent than 
heterogametic. A set of no-choice outcross experiments 
yielded only self-fertilized progeny from the Species A parent 
and reproductive failure from the pomilia parent, suggesting 
complete Species A x pomilia hybrid inviability. The third 
species of Physa inhabiting South Carolina, P. acuta , is more 
genetically similar to Species A but bears a distinctive penial 
anatomy. Mate choice tests uncovered no evidence of sexual 
isolation between Species A and P. acuta , and hybridization 
occurred readily, with some reduction in parental fecundity 
but normal FI viability. Species A x acuta FI hybrids appear, 
however, to be 100% sterile. Thus, the relationship between 
the degree of reproductive isolation and genetic divergence 
seems to be stronger than that between reproductive isolation 
and penial anatomy in the physid snails of South Carolina. 
Physa Species A warrants formal description. 

Additional keywords: Speciation, Physella , Physa acuta , Physa 
pomilia , mate choice, sexual isolation, hybridization, allozyme 
electrophoresis 


INTRODUCTION 

In recent years, a great deal of interest has foeused on the 
evolutionaiy biology of freshwater pulmonate snails in the 
family Physidae (Tsitrone et ah, 2003; Bousset et ah, 2004; 
I lenry et ah, 2005; 2006; Escobar et al. 2007). Their great 
reproductive plasticity, which includes selfing, mixed- 
mating, and outcrossing in either or both sexual roles, 
together with their ease of culture and the availability of 
genetic markers, lias made physid snails a favorite model 
for the study of sex allocation generally (Dillon and 


Wethington, 1992; Wethington and Dillon, 1991; 1993; 
1996: 1997). But despite great advances in our under¬ 
standing of broad aspects of their reproductive biology, 
progress in disentangling the complex evolutionary rela¬ 
tionships within the family Physidae has been slow. 

The classification system of George Te (1978; 1980) 
reeognized about 40 species and subspecies of physids in 
North America, arranged into genera and subgenera by 
penial anatomy. Within the group of nominal speeies 
bearing the penial eomplex Te characterized as “type-b,” 
however, Dillon and Wethington (2006a) reported no re¬ 
productive isolation among P. gyrina (Say, 1821), and five 
other more recently described speeies: P. ancillana (Say, 
1825), P. aurea (Lea, 1838), P. microstnata (Chamberlain 
and Berry, 1930), P. parken (Currier in DeCamp, 1881), 
and P. ntahensis (Clench, 1925). The addition of P. sayi 
(Tappan, 1838) to the list of type-b synonyms of P. gyrina 
was suggested by the survey of genetic variation at 
allozyme-eneoding loci offered by Dillon and Wethington 
(2006b). 

In the group of physids bearing Tes “penial eomplex 
type-c,” Dillon et al, (2002) could find no reproductive 
isolation among P Integra (Haldeman 1841) from the 
American northeast, P. lieterostropha (Say, 1817) from 
the American southeast, or the cosmopolitan P. acuta 
(Draparnaud, 1805), described from Europe prior to 
any American species of physid. Physa cubensis 
(Pfeiffer, 1839), from the Caribbean, and P. virgata 
(Gould, 1855), from the American West, have also re¬ 
cently been synonymized under P. acuta (Paraense N 
Pointier, 2003; Dillon et al., 2005). Reproductive isola¬ 
tion is complete, however, between physids bearing 
type-b and type-e penial complexes (Dillon et al., 2004). 

Te also reeognized a group with penial anatomy inter¬ 
mediate between type-b and type-c. These “type-be” 
species included P. henclersoni (Clench, 1925), originally 
described as a subspecies of P. pomilia (Conrad, 1834). 
But since Tes observations suggested to him that 
P. pomilia bore type-e penial morphology, he lowered 
pomilia to subspeeifie status under P. lieterostropha and 
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raised hendersoni to the rank of species. More thorough 
observations and experiments have confirmed, however, 
that topotypie P. pomilia bear penial anatomy type-bc, 
and that they are not reproductively isolated from 
P. hendersoni , a junior synonym (Dillon et ah, 2007). 

Recently a new classification has been proposed synthe¬ 
sizing laboratoiy experiments on reproduetive isolation 
together with mtDNA sequence divergence and morphol- 
ogieal observations (Wethington, 2004; Wethington and 
Lydeard, 2007). This classifieation recognizes approxi¬ 
mately 12 North American species and documents a loose 
correspondence between mtDNA sequence phylogroups 
and Tes penial morph ok )gies as outlined above. 

In addition, the sequence data of Wethington and 
Lydeard suggests that a previously unrecognized species 
of P/njsa , characterized by a dark body and elongated 
shell, might inhabit the swamps and ditches of the 
southeastern coastal plain. This species, bearing the 
type-be penial anatomy of P pomilia but genetically 
more similar to P. acuta , was referred to as “ Phijsa Spe¬ 
cies A”. The purpose of the present paper is to report 
the results of experiments designed to test for reproduc¬ 
tive isolation between Phijsa Species A and populations 
of the two other physids occurring in South Carolina, 
P. acuta (type-e) and P. pomilia (type-bc). 

The origin and evolution of reproductive isolation has 
been the subject of intense interest sinee the early twen- 
tieth-eentuiy birth of the Modern Synthesis (Mayr, 
1942; 1963). The barriers that may evolve between a 
pair of populations are conventionally divided into pre- 
zygotie components (such as sexual isolation) and post- 
zygotie components (such as hybrid inviability or 
sterility). The former is typically assessed using mate 
choice tests (Bateson. 1983) and the latter by no-ehoice 
breeding experiments (Coyne and Orr, 2004). Here we 
report the results of both mate choice and no-choice 
breeding experiments between a reference population 
of Phijsa Species A from South Carolina and P. acuta , 
then (separately) Phijsa Species A and P. pomilia. 


MATERIALS AND METHODS 

The Phijsa acuta population used to found “line A” for 
these experiments inhabits the main pond at Charles 
Towne Landing State Park, west of the Ashley River, 
within the city limits of Charleston, SC (32.8062° N, 
79,9862° W). Snails of this population are not reproduc¬ 
tively isolated from P. acuta sampled near the type local¬ 
ity for the species in France (Dillon et al., 2002). The 
Phijsa pomilia population used here to found “line IF 
was collected from the type locality for Phijsa pomilia 
hendersoni (Clench. 1925): the Combahee River at the 
US 21/17A bridge, 1 km E of Yemassee, Hampton 
County, SC (32.7060° N, 80.8281° W). Dillon et al. 
(2007) reported no reproductive isolation between this 
population and snails sampled from Conrads (1834) 
type locality for Phijsa pomilia sensu stricto in Alabama. 
The reference population of Phijsa Species A used to 


found line “S” was collected from the spring by Huger 
Creek at Huger Landing, 4 km N of Huger, Berkeley 
County, South Carolina (33.1305°N; 79.8111°W). 

All snails were cultured in transparent polyethylene 
10 ounce drinking cups filled with approximately 
210 ml of aerated, filtered pond water and covered with 
a 95 x 15 mm polystyrene Petri dish lid. They were fed 
O.S.I. Spirulina Aquarium Flake Food, sold in pet stores 
primarily as a diet for herbivorous aquarium fishes. All 
experiments took place at room temperature, approxi¬ 
mately 23°C. I initially isolated ten wild-collected snails 
from each study population in separate cups, collected 
egg masses with weekly water change, and reared the 
offspring to 2 mm shell length, approximately 3 weeks 
post-hatching (well in advance of maturity). These three 
sets of wild-collected but laboratory born sibships were 
designated Al through A10, SI through S10, and III 
through 1110. From these sibships were drawn isolates 
for the mate choice tests and pairs of parents for the 
study of postzygotie reproduetive isolation. 

For mate choiee tests, large samples of juvenile snails 
from all three populations were reared to maturity 7 over 
the course of 8—10 weeks isolated in individual cups, 
with weekly feeding and water change. Two experiments 
were performed: one comparing Speeies A to P. acuta 
and the other comparing Speeies A to P. pomilia. Each 
experiment was composed of three trials, eaeh trial in¬ 
volving 10 adult snails from one population and ten adult 
snails from a second, all approximately matching in their 
shell sizes. Snails were blotted dry and marked with a 
small dab of fingernail polish according to their popula¬ 
tion of origin. Then the 20 individuals were simulta¬ 
neously introduced into a 2 liter glass beaker (filled 
with 1,400 ml ol filtered, aerated pond water) and 
placed on a glass table to facilitate observation. 

Mating activity was monitored for 6 hours. When a 
snail first successfully copulated as male (defined as the 
complete insertion of its penis into the gonopore of a 
partner) it was removed from the beaker, its shell marked 
with a dot of white correction fluid, and returned. Eaeh 
individual was often involved in many matings over the 
6 hours of observation, both in the male and in the female 
role, but only its first successful copulation in the male 
role was recorded. This was an arbitrary decision on my 
part (since both eopulants in a pair might mate in either 
role, and the result is not a “choice” but rather the out¬ 
come of a contest), but neeessary nevertheless to prevent 
double-counting. Note that this design yields a slight bias 
toward heterogametic pairings, not 1.1 but rather 9:10. 

Each trial involved 20 fresh snails, entirely unmated. 
Three such trials were performed testing for sexual iso¬ 
lation between Species A and acuta (the SA experiment) 
and three additional trials performed testing for sexual 
isolation between Speeies A and pomilia (the SH experi¬ 
ment), pooling results within experiment to yield a maxi¬ 
mum of 60 observations in each ease. Chi-square 
statistics were calculated from the pair of 2 x 2 contin¬ 
gency tables that resulted, normalized by 4/N, as a mea¬ 
sure of sexual isolation (Gilbert and Starmer, 1985). 
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For no-choice tests of postzygotic reproductive isola¬ 
tion, three sets of incross control cups were established 
using pairs ol unrelated parents drawn from the ten sib- 
ships within each of the populations (S, H, and A), as for 
example SlxS2, S2xS3, ..., SlOxSl. Two sets of out- 
cross experimental cups were also established with 10 
pairs of snails across populations, the SA cross (SxAl, 
SxA2, ..., SxA10) and the SH cross (SxIII, SxH2, ... 
SxIIlO). Each pair of parents received a water change 
and fresh food eveiy 7 days, at which time the sides of 
the cup were inspected for egg masses. (Note that any 
egg mass might result from outcrossing, or be the prod¬ 
uct of self-fertilization by either parent.) If egg masses 
were present, all embryos were counted and adults 
transferred to a fresh cup. Eggs were monitored until 
hatching (generally about 2 weeks) and all viable, crawl¬ 
ing FI juveniles counted. Observation was terminated 
upon the death of either parent in a pair. 

Crosses were initiated with pairs of snails aged one 
week post hatch. Then any difference in the central 
tendency of age at first reproduction (in weeks post 
hatch) between the 10 outcross pairs and the combina¬ 
tion of both sets of 10 corresponding control pairs was 
tested by calculating a combined (30 pair) median and 
comparing counts above and below that median using 
Fishers exact tests. 

For statistical analysis of fecundity and FI viability, 
week 1 was established separately for each set of 10 
pairs as the first week in which eggs were laid by 3 or 
more pairs of parents. Embryos and viable hatchlings 
were subsequently counted for 10 weeks. I then aver¬ 
aged the embryo production of each pair of parents 
across its lifetime, ignoring any leading (pre-maturity) 
zeros and any postmortem zeros, while including as 
0 any failure to reproduce by viable, mature pairs. So, 
for example, if one parent in a pair of snails died at week 
6, leaving a record of 0, 0, 40, 0, 50 embryos for the pair, 
their mean fecundity would be 90/3 = 30 embryos per 
week. A Kruskal-Wallis nonparametric ANOVA was used 
to test whether any significant difference existed in the 
central tendency of weekly mean fecundity of either set 
of 10 outcross pairs (SH or SA) and the 2 corresponding 
sets of 10 control pairs. 

Similarly, I averaged the counts of FI hatchlings with¬ 
in pairs across weeks, ignoring zeros not corresponding 
to embryo production, and divided each pair mean by its 
mean embryo production to obtain pair mean FI viabili¬ 
ty. If 35 + 45 hatchlings were recovered from the exam¬ 
ple pair of snails above, their mean FI viability would be 
(35/40 + 45/50)/2 = 88.9%. A second Kruskal-Wallis non¬ 
parametric ANOVA was used to test whether any signifi¬ 
cant difference existed in the central tendency of weekly 
mean FI viability posted by either set of 10 outcross pairs 
and its 2 corresponding sets of 10 control pairs. 

To assess the fertility of putative hybrid offspring, FI 
hatchlings (from both experimental sets and all three 
control sets) were reared from each of 3 separate unre¬ 
lated pairs to size 2 mm. These were crossed in time 
series: 1 early pair from eggs laid around week 1, 1 mid¬ 


dle pair produced around week 5, and 1 late pair pro¬ 
duced around week 10, to yield 9 FI pairs. So if the 
putative hybrid progeny were reared from pairs SxAl, 
SxA2, and SxA3, for example, they were crossed as 
SAlxSA2 early, SA2xSA3 early, SA3xSAl early, 
SA1 xSA2 middle, SA2xSA3 middle, ..., SA3xSAl late. 
Nine crosses were likewise constituted for corresponding 
controls S and A, and the total of 3 x 9 = 27 crosses of 
FI snails reared to adulthood for each experiment, with 
weekly feeding and water change. An identical set of 
27 cups was established to evaluate hybrid fertility in the 
SH experiment. I recorded the dates at which embryos 
and viable F2 hatchlings were produced by each pair. 

A larger sample of FI progeny from 3 outcross pairs 
from both the SA and SH experiments were reared to 
4-5 mm shell length, at which time they were frozen in 
100 pi of tissue buffer for analysis by allozyme electro¬ 
phoresis. We have identified 12 enzyme-encoding loci at 
which allozyme variation is interpretable as the product 
of codominant alleles segregating in Mendelian fashion 
(Dillon and Wethington, 1994). These are aconitase 
(Aeon), esterases (three loci: Estl, Est3, Est6), glucose 
phosphate isomerase (Gpi), isocitrate dehydrogenase 
(two loci: Isdhl and Isdh2), leucine aminopeptidase 
(Lap), mannose phosphate isomerase (Mpi), phospho- 
glucomutase (two loci: Pgml and Pgm2), and 6-phospho- 
gluconate dehydrogenase (6pgd). We used horizontal 
starch gel electrophoresis in an aminopropylmorpholine 
pH 6 buffer system to resolve allozyme variation at the 
Gpi, Tsdli, and 6pgd loci, a Tris-Citrate pH6 buffer sys¬ 
tem for Aeon, Mpi, and Pgm, and a TEB8 system for 
6pgd, Lap, and Est. Details regarding our electropho¬ 
retic methods, including a description of our equipment 
and recipes for stains and buffers, have been previously 
published (Dillon, 1992; Dillon and Wethington, 1995). 

The set of no-choice mating experiments described 
above were conducted simultaneously with those of Dil¬ 
lon et al. (2007), using identical techniques. The data 
reported here on the reproductive performance of the 
A and H incross control lines have been published pre¬ 
viously, although the SA and SII experimental results, as 
well as the S incross control, are original to the present 
investigation. 


RESULTS 

The SA mate choice experiments did not reveal any 
evidence of sexual isolation between Species A and 
P. acuta (Tabic 1, upper). A total of 49 copulations were 
observed (of a possible 60 total), apparently without 
regard to species (normalized x^ = 0.82, p = 0.37). The 
SI I experiments did, however, suggest prezygotic repro¬ 
ductive isolation between Species A and P. pomilia 
(Table 1, lower). The 38 copulations observed in the SH 
mate choice tests included only 2 of pomilia inseminated 
by a Species A partner, while 10 pomilia were insemi¬ 
nated by pomilia partners. There was also a bias toward 
homozygotic pairings on the Species A side, yielding a 
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Table 1 . Copulations observed in the two mate choice 
experiments, Plujsa Species A x P. acuta (above) and Plujsa 
Species A x P. pointful (below). 



Males 


Homogametic 

Ileterogametic Totals 

Females Species 

15 

15 30 

A (S) 

P. acuta (A) 

12 

7 19 

49 

Females Species 

16 

10 26 

A (S) 

P. pomilia (H) 

10 

2 12 

38 


significant overall deviation from random mating (nor¬ 
malized x 2 = 6.63, p = 0.01). 

Reared together in a no-choice design, mixed pairs of 
Species A and P. acuta showed no delay in age at first 
reproduction, their modal age at maturation (7 wks) 
indeed slightly less than that observed in either matched 
Species A or matched acuta control pairs (Table 2). 
A reduction was apparent in parental fecundity, however, 
tilt 1 median of 55.2 embryos/wk posted by SA outcross 
pairs significantly below both controls (p = 0.027). 
The 73.1% median viability of the FI Species A / acuta 
hybrids was intermediate between the FI viabilities 
observed from incross controls. 

Electrophoretic analysis of a sample of offspring from 
three SA outcrosses confirmed the hybridity of all FI 
progeny. One pair of parents was fortuitously fixed for 
alternative alleles at the Isdh locus, yielding a sample ol 
twelve entirely heterozygous progeny. A second pair of 
SA parents were both heterozygous at the Est3 locus 
(Est3 100 /Est3 96 x Est3 96 /Est3 92 ), yielding twelve FI proge- 
ny in four classes. The third pair of parents included one 


Table 2. Statistics comparing the fitness of Plujsa Species 
A x P. acuta outcrosses to pure Phifsa Species A and pure 
P acuta controls. 



Species A 

SA outcross 

P. acuta 

First opposition, P generation 

Mode (weeks post hatch) 8 

7 

8 

Range 

7-8 

7-10 

5-10 

Weekly mean parental 
fecundity 

Median (embiyos) 

75.4 

55.2 

66.9 

Range 

19-105 

18-81 

17-104 

Weekly mean FI viability 

Median (%) 

80.5 

73.1 

61.7 

Range 

66-91 

49-99 

45-86 

FI Fertility 

67% 

0% 

100% 

Viable F2 hatch 

Median (week) 

12.5 

- 

9 

Range 

4-19 

- 

9-11 


heterozygote at the Isdh locus (Isdh l()0 /Isdh 1(K) x Isdh 100 / 
Isdli 07 ), and one heterozygote at the Est3 locus (Est3 100 / 
Est3 92 x Est3 92 /Est3 92 ), yielding at both loci twelve FI 
progeny representing the heterozygous and one homozy¬ 
gous class, missing the other homozygous class entirely. The 
likelihood of missing a single homozygous class in twelve 
selfed progeny from a heterozygous parent would be 0.032. 

None of the nine pairs of FI progeny from the SA 
outcross produced viable F2 offspring. One SA pair was 
terminated early by mortality, while the other eight pairs 
all laid eggs profusely, beginning at week 7 and extend¬ 
ing to week 19. All egg masses laid by all eight pairs of 
SA hybrids over the 12 week period were held for five 
weeks, with no hatching observed. 

Reared together in a no-choice design, pairs of Species 
A and P. pom ilia demonstrated significant delays in age at 
first reproduction behind that posted by their combined 
controls (Fishers exact p = 0.003). Their modal age ol 9 
weeks at the onset of egg laying was slightly behind both 
the Species A control and the pomilia control (Table 3). 
The median parental fecundity of 27.3 embryos/wk 
posted in the SI I outcross experiment was also signifi¬ 
cantly lower than both incross controls (p = 0.002), and 
the median viability ol their progeny (64.S%) lower than 
the Species A control. Only one of the nine pairs of first 
generation progeny from the SII experiment yielded via¬ 
ble second generation offspring, at week 20. Most of the 
remaining first generation pairs laid eggs that failed to 
hatch, generally over many weeks ol observation. 

Electrophoretic analysis revealed that two sets of Sll 
parents were fortuitously fixed for alternative alleles at 
the LAP locus. Samples of ten first generation progeny 
from both of these crosses yielded only one homozygous 
class, strongly suggesting self-fertilization by the Species 
A parent, and no reproduction by the pomilia parent. 
Absence of suitable genetic markers made inference re¬ 
garding the third set of SII progeny analyzed equivocal. 


Table 3. Statistics comparing the fitness of Plujsa Species 
A x P. pomilia outcrosses to pure Plujsa Species A and pure 
P. pomilia controls. 



Species A 

SH outcross 

P. pomilia 

First opposition, P generation 

Mode (weeks post hatch) 8 

9 

7 

Range 

7-8 

8-12 

7-11 

Weekly mean parental 
fecundity 

M edian (e mb nos) 

75.4 

27.3 

57.4 

Range 

19-105 

7-41 

19-68 

Weekly mean FI viability 

Median (%) 

80.5 

64.8 

55.8 

Range 

66-91 

10-86 

33-86 

FI Fertility 

67% 

10% 

78% 

Viable F2 hatch 

Median (week) 

12.5 

20 

5 

Range 

4-19 

- 

3-8 
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DISCUSSION 

The experiments reported here confirm reproductive 
isolation between the “ Plujsa Species A” of Wethington 
(2004) and populations representing both of the other 
physid species inhabiting South Carolina, P. acuta and 
P. pom ilia. An initiative to formally describe Species A 
has just been published (Wethington et al., 2009). The 
reproductive isolation displayed by these three species is 
of different degrees, however, and apparently more 
closely related to their genetic divergence than to their 
reproductive anatomy 

Plujsa Species A and P. acuta cluster in the same 
mtDNA phylogroup (Wethington and Lydeard, 2007) 
but differ in their penial anatomy. The mate choice tests 
reported here yielded no evidence of sexual isolation 
between them. A significant reduction in the joint fecun¬ 
dity of Species A x P. acuta outcross pairs was indeed 
revealed by no-choice breeding experiments, although 
there was no evidence of reduced viability in the FI 
hybrids such crosses produced. Species A x acuta 
hybrids were, however, entirely sterile. 

Plujsa Species A and P. pomilia share identical type-bc 
penial anatomy, but are more distantly related genetical¬ 
ly. The reproductive isolation that Species A and pomilia 
display under controlled conditions is of a greater degree 
than that observed between Species A and acuta. Paired 
in a no-choice design, Species A and pomilia parents 
displayed delayed reproduction, reduced fertility, and re¬ 
duced offspring viability. All the viable first-generation 
progeny recovered from SH outcrosses were attributable 
to self-fertilization by the Species A parent, reproduction 
by the pomilia parent apparently foreclosed. It seems 
likely that the substantial reduction in reproductive suc¬ 
cess demonstrated in the second generation by SH off¬ 
spring may be attributable to Species A inbreeding 
depression. Such results are quite simiar to those we 
obtained from crosses between the type-c Plujsa acuta 
and the type-b Plujsa gijrina (Dillon et ah, 2004). 

In addition, mate choice tests returned evidence of 
prezygotic reproductive isolation between Species A 
and P. pomilia. Physids seem to mate according to a 
modified “Batemans Principle” (Bateman, 1949). They 
arc generally quick to copulate as males, and display 
little discrimination, but when mounted in the female 
role they can be choosy, often displaying rejective be¬ 
haviors like evasion and shell-shaking (DeWitt, 1991: 
1996; Wethington & Dillon, 1996; McCarthy and Sih, 
2008). Only 38 of the 60 snails tested in the SI I mate 
choice experiments were ultimately able to mate as 
males, at least partly because of the high frequency of 
rejective behaviors they encountered in heterogametic 
couplings. Our observation that pomilia copulants were 
more rejective of heterogametic insemination than Spe¬ 
cies A copulants may be related to our separate observa¬ 
tion that the Species A partner in our SH no-choice 
experiments retained the ability to reproduce by self- 
fertilization, while the pomilia partner apparently did 
not. This situation is similar to that we have previously 


described for the interspecific pairing of Plujsa acuta 
and P. pomilia (Dillon et al., 2007). 

Although nothing is known about the genetics of re¬ 
productive isolation in pulmonate snails, a large body of 
research has confirmed that both prezygotic and postzy- 
gotic barriers are inherited in a complex and polygenic 
fashion in Drosophila (Wu and Palopoli, 1993; Ritchie 
and Phillips, 1998). In general it has been found that 
postzygotic isolating mechanisms evolve independently 
of, and tend to lag behind, prezygotic mechanisms 
(Coyne and Orr, 2004). Whether the latter can be rein¬ 
forced by natural selection on the former is controver¬ 
sial. Experiments to trace the evolution of both sets of 
characters through the larger phylogeny of the Physidae 
are currently ongoing. 
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